Detecting Entanglement with Jarzynski's Equality 
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We present a method for detecting the entanglement of a state using non-equilibrium processes. 
A comparison of relative entropies allows us to construct an entanglement witness. The relative 
entropy can further be related to the quantum Jarzynski equality, allowing non-equilibrium work to 
be used in entanglement detection. To exemplify our results, we consider two different spin chains. 
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In quantum information theory, entanglement is con- 
sidered not only an interesting phenomenon, but also a 
resource which can be used in quantum computation. 
Entanglement has therefore been the topic of much re- 
search. A separable state can be written as a convex 
sum of pure product states, a = pi<jj of • • • erf 
where the piS are the weights of the product states, of, 
with J2 i Pi = 1, while an entangled state cannot. Many 
methods have been devised to measure and detect en- 
tanglement, even for thermal and for many-body sys- 
tems [![. The entanglement witness [2] is an expectation 
value of an operator which is bounded for any separable 
state, whereas entangled states can exceed this bound. A 
thermodynamic witness allows us to use thermodynamic 
quantities such as the magnetic susceptibility [3] to de- 
tect entanglement. The major advantage of using a such 
a witness is that we can detect thermal many-body en- 
tanglement using experimentally measurable quantities. 

Thus far, these thermodynamic witnesses have only 
been used for detecting entanglement in equilibrium sys- 
tems. However, a result from condensed matter theory, 
Jarzynski's equality 0], allows the change in free energy 
between two equilibrium states to be related to the non- 
equilibrium work done needed to drive the system from 
one state to the other. While the work done can be mea- 
sured or calculated in an experiment, the change in free 
energy cannot. Thus Jarzynski's equality can be used to 
experimentally estimate the change in free energy dur- 
ing a non-equilibrium process fB|. It is the aim of this 
letter to use Jarzynski's equality to witness equilibrium 
entanglement using non-equilibrium processes. Our work 
also raises the exciting possibility of using this witness to 
detect entanglement in biological systems. 

In our construction of an entanglement witness, we use 
the relative entropy, a directed distance from an initial 
state pi to a final state p/, given by 



S (pf\\pi) = tr (Pf lo gP/) ~ tr (pf log Pi). 



(1) 



for the relative entropy of entanglement, Ere(p/) — 
mm Pie s S(pf\\pi), where we take the minimum over the 
set of separable states <S, to find a css . The relative en- 
tropy can measure entanglement for both equilibrium 
and non-equilibrium, pure and mixed states, and there- 
fore for thermal, open and closed systems. 

We can now construct an entanglement witness using 
the relative entropy by introducing an arbitrary state p* . 
Since the set of separable states is convex, and the rela- 
tive entropy is a directed distance, if the distance from 
<j css to p is larger than the distance from p* to p, then 
p* is entangled. Fig. ([T]) gives a two dimensional repre- 
sentation of this idea. Hence our witness is 



S( P \\a css ) > S(p\\p*). 



(2) 



If p* satisfies this inequality, we know it must be entan- 
gled. The witness is best when p is a pure state and hence 
is located at the edge of the outer ellipse in Fig. ([T]). We 
will refer to this inequality as the relative entropy wit- 
ness. We note that this witness can detect entanglement 
in both equilibrium and non-equilibrium states. 

Although originally a classical result, it has been shown 
that Jarzynski's equality, 



/3W\ e -/3AF 



(3) 



This is a measure of entanglement [6( when pi = a css 
is the closest separable state to pf. Formally, we have 



where f3~ l is the temperature, W is the work done on the 
system and AF is the change in free energy between the 
initial and final equilibrium states, is valid for both open 
Q and closed quantum systems [I, 0, [To[ • The brackets 
(• • •) denote an average over all possible realisations of 
the work, or trajectories in phase space. Both the path 
and the rate at which the system is driven are fixed for 
the equality, though each are arbitrary. 

There are several different methods (for a review, see 
reference [TT| ) , in the literature for deriving a quantum 
version of Jarzynski's equality, however we discuss the 
one which has been successfully theoretically verified [13] • 
In a closed quantum system, instead of classical trajec- 
tories in phase space, we define the quantum equivalent 
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FIG. 1: This is a 2d representation of the multidimensional 
Hilbert space. The small oval represents the set of separable 
states, and the large oval the set of all states. a cm is a com- 
pletely mixed state, and a css is the closest separable state to 
p. Any state along the pink curve, such as a' , has the same 
"distance" in terms of the relative entropy as (j ess to p. 



of quantum transition probabilities [1, 0, [Io[. An ini- 
tial Hamiltonian Hi and a final Hamiltonian Hf have 
eigenvalues E % n1 E^ and eigenvectors |<^), respec- 
tively. We perform a measurement of the energy at time 
U and then again at tf so that the system is in a spe- 
cific energy eigenstate. The quantum transition prob- 
abilities are then defined as q myn = |(^4l^(^/)l^n)| 2 
where U(tf) = T < e~*^o / H ( s ) ds [ s the time evolution 
operator, and T< is the time ordering operator. q m ,n 
can be interpreted as the probability that the final state 
of the system is |(/>4) given that it was initially in the 
state \(j) l n ). The average is then given as (e _/3W ) = 
En( e ~ PK /^)ErnQrn,ne-^ E -- E ^ where the work is 
defined as W = E^ - E l n and Z { = tr{e~^ Hi ) is the 
initial partition function. 

Consider now an open quantum system (subsystem, 
S) interacting with a bath, £?, with total Hamiltonian 
H(t) = Hs(t) + Hsb + Hb and arbitrary coupling, Hsb 
0]- As only the subsystem is time dependent, the change 
in energy of the total system equals the work done on the 
subsystem. Hence the average in equation (j3]) is identical 
to the closed system case. Further, the free energy of the 
total system is given by F(t) = Fs(t) + Fb- This allows 



Jarzynski's equality to be written (e " w ) 



-PAF S ^ 

Other fluctuation theorems have also been derived. 
One equality which will be useful Q is ( e -(PfEf-PiEi)} = 
e -(Pf F f-Pi F i) m This demonstrates that a change in tem- 
perature between the initial and final state can also be 
taken into account. However, unless = /?/, the quan- 
tity (PfEf — PiEi) no longer relates to work. We refer to 
this equation as the Jarzynski-Tasaki equality. 

Consider again the relative entropy, equation (pQ). We 
now restrict each of the states to be in thermal equi- 
librium. Thus we have initial and final states, a css = 
e -PiHi/Zi and p = e~@ fHf /Zf respectively. Similarly, 
p* = H /Z*. Expanding equation (pQ), we can write 
the relative entropy in terms of a change in free en- 
ergy [3, [3, S(p\\a css ) = A(pF) - tr(p f A(f3H)) where 



A(pF) = /3fFf-/3iFi and A(pH) = P f H f -piHi. Com- 
bining this identity with the Jarzynski-Tasaki equality, 
we find that 

S(p\\a css ) = -tr( P A({3H)) - ln< e -^-^>). (4) 

Equation (|4]) relates the entanglement to the average 
change in energy at different temperatures (in a possi- 
bly driven system). When Pf = we can instead relate 
the entanglement to the average work done in creating 
the quantum correlations of p from the purely classical 
correlations of cr css . In addition to this being an inter- 
esting result in itself, we can also use this definition of 
the relative entropy in the entanglement witness, equa- 
tion ([2j). We call this the relative Jarzynski witness and 
discuss this in more detail below. 

An open quantum system (or subsystem) can also 
be considered. As discussed previously, such systems 
also obey Jarzynski's equality, (e _/3W ) = e~@ AFs where 
Fs(t) is the free energy of the subsystem. We define 
Y(t) = tr{e-^ Hs{ - t)JrHsBJrHB) ) as the partition function 
of the total system and Zb = tr(e~^ HB ) as the partition 
function of the bath. The partition function of the sub- 
system, Zs(t) = tr(e~P Fs ) = Y(t)/ZB can be associated 
with an effective Hamiltonian [7[, 



H*ff(t) = --In 



tr B (e 



-P(H s (t)+H S B+H B ) 



tr B (e-P HB ) 



(5) 



so that Zs(t) = trs(e ^ Re ^). Using these equations, 
and since the initial and final states must be in equi- 
librium, we have ps = e~ f3Heff ^/Zs(t) for each state. 
This only represents the state of the system when it is in 
equilibrium. 

The relative entropy can now be defined in terms of the 
effective Hamiltonian and the work done on the subsys- 
tem, S(ps\\as, css ) = -f3tr[p s (H e / f - Hf f )] -\n(e^ w ). 
Hence we can write the entanglement witness for an 
open quantum system, S(ps\\o's,css 

) > S(ps\\p* s ), where 



Ps 
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/Zg as the state is in equilibrium. 



We have shown that it is possible to detect entangle- 
ment in a state p* or p* s using a non-equilibrium process. 
There are three ways we can use the entanglement wit- 
ness. First, if we have a specific state p* in mind but 
don't know whether it is entangled, the witness allows us 
to detect entanglement in this state. Second, if we have 
the Hamiltonian of a system, we can detect entanglement 
in that system. We ask for which values of the parame- 
ters of the system, such as a magnetic field, is the system 
entangled? In this case we find out in which p*s of the 
system entanglement can be detected. Computationally, 
linking equation ([2]) to Jarzynski's equality simply gives 
a different way to calculate the relative entropy. It is 
in the third method, the experimental applications, that 
are exciting in this respect. 



Experimentally, we can relate the relative Jarzynski 
witness to non-equilibrium processes. For the closed 
quantum system when f3f = and the open quantum 
system, this corresponds to a series of measurements. We 
drop the subscript S that denotes the subsystem in the 
open quantum system here since the work done in the 
open system is equal to the change in energy of the to- 
tal, closed, system. Hence the discussion is valid for both 
open and closed systems. 

Since we consider a quantum system, measurements 
of the energy on many replicas of the same system will 
give different values. As the quantum Jarzynski equality 
demonstrates, each time we measure an initial and a final 
energy of a system to calculate the work, we obtain differ- 
ent results. After many measurements of the initial state 
(a C ss) and the final state (p), we can calculate the average 
(e _/3W ). We then repeat this procedure with initial state 
p* and compare the resulting experimental values of the 
relative entropy. The value of tr(p(Hf — Hi)) can also 
be experimentally measured. For instance, if a magnetic 
field is driving the process, this corresponds to the change 
in the field multiplied by the final state magnetisation. 
We can now detect entanglement in p*. 

When f3f ^ we can use the Jarzynski- Tasaki equal- 
ity and a similar argument holds. However, it is no longer 
the work done that is measured. Instead we measure the 
initial and final temperatures of the system in addition 
to the energy eigenvalues. 

A problem with using our Jarzynski witness is the pos- 
sibility that a css is not an equilibrium state of the sys- 
tem, and therefore we cannot define Jarzynski's equal- 
ity. However, we find that we do not require <j css to 
be in equilibrium itself. Instead, we require only that 
we have an equilibrium state a' of the system where 
S(p\\(T css ) = S(p\\cr f ). The states satisfying this equal- 
ity are represented by the pink curve in Fig. (pQ). 

We now illustrate the entanglement witness with two 
examples. We first consider a three qubit XX Z spin 
chain as we can define both the initial and final states 
to be in equilibrium. In the second example, we use a 
seven qubit chain to demonstrate what happens when 
the closest separable state is not in equilibrium and we 
must use a different equilibrium state a' . 

For each example, we have calculated the witness us- 
ing the relative entropy witness and using the Jarzynski- 
Tasaki witness, and find both give the correct results. 
This also allows us to successfully numerically verify 
the Jarzynski-Taskaki equality. We calculate the time 
evolution operator exactly in the three qubit case as 
\H(ti), H(t2)\ = 0, and using the method described in 
[151 ] for seven qubits as [H (ti) , H fe)] ^ 0. This method 
allows an approximation of U(tf) to be calculated us- 
ing U(t f ) = Un=o le ~ iH{tn)At - We use At = °- 001 to 
give accurate results. We use these examples rather than 
that of an open quantum system since the closest sepa- 
rable state to p n given below is known. This allows us 




FIG. 2: This plots J z versus B and T when N = 3, and 
shows the values for which we can detect entanglement in p* . 
The state is entangled between the surface of the plot and the 



to do some of the calculation analytically which allows 
further insight into the problem. 

We take our state to be close to the pure symmetric 
state, p n = (l/n)5(|00---01))5((00---01|) where S is 
the total symmetrisation operator, whose closest separa- 
ble state [HI, [TtJ is known to be 



n — / (n~ 



i) fc 5(|e...m))S((e...m|). 

k n — k k n — k 

(6) 

We will identify the states p n and a C ss,n with thermal 
equilibrium states, e~@ H jZ, and hence we will not have 
exactly the states above. However, we find that the rela- 
tive entropy calculated in each case is identical to many 
significant figures. 

The Hamiltonian of the XX Z spin chain is 



H 



E 



2 \ a i 



l+l 



■Baf 



(7) 

where J and J z are coupling strengths, and B is a mag- 
netic field. 

For our first example, the three qubit spin chain, 
we require J z and B to be time dependent. Both ps 
and a C ss,3 can be written as thermal equilibrium states, 
p = e~@ H jZ, of the Hamiltonian. For the initial state 
to be <t CSS5 3, we require that B\ = /3 _1 log(2)/2 and 
J Zi i = (2 J — log(3))/4 at a low temperature. For 
concreteness, we take = 0.01 and J = 1. For the 
final state to be p3, we have Bf = 1/2 and J z j = 0. 

We can now detect entanglement in an arbitrary equi- 
librium state, P3 using the entanglement witness. Fig. 




A possible application of this work is the detection of 
entanglement in biological systems. The photosynthetic 
bacteria, Prosthecochloris aestuarii, can be modelled us- 
ing a seven spin Hamiltonian. Using experimental values 



FIG. 3: This plots J z versus B and T when N = 7, and 
shows the values for which we can detect entanglement in p* . 
The state is entangled between the surface of the plot and the 



© shows the values of the magnetic field, J z and the 
temperature for which we can detect entanglement: we 
can detect that p^ is entangled in the region between 
the surface and the axes. Hence, experimentally driving 
a system from the state p% with values of 5, J z and T 
that are within the surface to the state p will allow en- 
tanglement to be detected on comparison with the same 
process starting at <t CSS) 3. 

Our second example is the 7 qubit spin chain, with 
B time dependent and J z = 0. For any 7 qubit 
chain, we cannot identify cr css j with a thermal equi- 
librium state, and hence we use <j f 7 = [(7 7 — 7 x 
6 6 )|0000000)(0000000|+6 6 6'(|0000001))6'((0000001|)]/7 7 
instead. For the initial state to be cr 7 , we require that 
Bi = p- 1 log[70993/46656]/2 + J at a low temperature. 
For concreteness, we again take f3~ l = 0.01 and J = 1, 
and for the final state to be p 7l we have Bf = 0.92. 

We can now detect the entanglement of a state p 7 as 
before. Fig. ([3]) shows the values of J z and the tem- 
perature for which we can detect entanglement. Again, 
we can detect that p 7 is entangled in the region between 
the surface and the axes. 

We note that although J z = for both the initial and 
final state Hamiltonians, this is not necessarily so for p*. 
Indeed, we can detect when p* is entangled in many other 
situations. This is due to the fact that the Hilbert space 
of the Hamiltonian is spanned by the set of n compu- 
tational eigenvectors, {|00 • • • 0), |00 • • • 01) • • • |11 • • • 1)}. 
Hence the entanglement witness applies to any state p* 
that exists within this Hilbert space. For example, we 
could introduce a Dzyaloshinskii-Moriya interaction to 
the Hamiltonian of p* and still use the witness to detect 
entanglement in the system. 



18l ll9[ and simplifying the model to an isolated system, 
we can use this Hamiltonian to construct a specific state 
p*. The 7 qubit chain defined above can then be used 
in the witness. In this simplified model, we do not de- 
tect any entanglement. However, we expect that the full 
model, and a more appropriate Hamiltonian Hf which is 
closer to H* will allow entanglement to be detected. 

We have presented a witness which uses the relative 
entropy to detect entanglement. When the states are in 
equilibrium, we have shown that Jarzynski's equality can 
be used to detect entanglement. Hence this witness en- 
ables entanglement to be detected using non-equilibrium 
processes. Using this witness, we have considered two 
examples. In one we can define an equilibrium closest 
separable state to p, and in the other we instead define 
an entangled equilibrium state which has the same di- 
rected distance to p in terms of the relative entropy. 
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